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NASA Langley Research Center, Hampton, Virginia 
Peter G. Iiju ; 

University of Florida, Gainesville, Florida 

Abstract 

Micro aerial vehicles have been the subject of continued interest and 
development over the last several years. The majority of current vehicle 
concepts rely on rigid fixed wings or rotors. An alternate design based on an 
aeroelastic membrane wing has also been developed that exhibits desired 
characteristics in flight test demonstrations, competition, and in prior 
aerodynamics studies. This paper presents a simulation model and an 
assessment of flight control characteristics of the vehicle. Linear state space 
models of the vehicle associated with typical trimmed level flight conditions 
and which are suitable for control system design are presented as well. The 
simulation is used as the basis for the design of a measurement based nonlinear 
dynamic inversion control system and outer loop guidance system. The 
vehicle/controller system is the subject of ongoing investigations of 
autonomous and collaborative control schemes. The results indicate that the 
design represents a good basis for further development of the micro aerial 
vehicle for autonomous and collaborative controls research. 


Introduction 

Micro aerial vehicles, or "MAVs'’, are typically 
designated as a class of aircraft with a maximum 
dimension of 6 inches that are capable of operating at 
speeds of 25 mph or less. 1 11 Developments in 
miniaturized digital electronics, communications, and 
computer technologies and strong support by DARPA 
have moved the prospect of very small autonomous 
flight vehicles from the realm of science fiction to 
science fact. The goal is for these vehicles to provide 
inexpensive and expendable platforms for surveillance 
and data collection in situations where larger vehicles 
are not practical. For example, they can be used for 
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battlefield surveillance or mapping the extent of 
chemical/radiation spills or viral outbreaks. Other 
applications include use in search and rescue 
operations, traffic/news coverage, and crop or wildlife 
monitoring. Many potential uses would require 
cooperative and collaborative control capabilities so 
that large numbers of MAVs could be used to cover a 
large operational area. In these types of applications 
MAVs could be coordinated from a central base 
station or used in collaborative swarms to collect and 
transmit data. 

The research and development required for 
developing MAVs and related systems is technically 
challenging and requires a number of technological 
advances that may benefit a broad range of aerospace 
applications. The development of a vehicle could also 
foster development of component technologies and 
help to support an emerging growth market for micro 
aerial vehicles. 

An aeroelastic fixed wing micro aerial vehicle 
concept has been developed by a team at the University 
of Florida with a goal to design a vehicle that could 
win the ISSMO (International Society of Structural 
and Multidisciplinary Optimization) Micro Aerial 
Vehicle Competition; a goal that has been 
accomplished each of the last four years. 1231 

The vehicle exploits an innovative aeroelastic wing 
with the ability to adapt to atmospheric disturbances 
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Figure 1 - photograph of Univ. of Florida MAV. 


and provide smoother flight thereby providing a better 
surveillance platfonn and making the vehicle easier to 
fly. This is accomplished via a passive adaptive 
washout mechanism. 

The adaptive washout technique has been taken from 
sailing vessels which use sail twist to greatly extends 
the wind range of the sail and produce more constant 
thrust (lift) in gusty wind conditions. Adaptive 
washout is produced in the MAV by deformation of the 
membrane wing in response to changes in speed and 
vehicle attitude. The result produces changes in wing 
camber and angle of attack along the span. The effect is 
to reduce the sensitivity of the vehicle to disturbances. 

NASA is collaborating with the University of 
Florida to develop an understanding of the underlying 
physical phenomena associated with the vehicle 
concept with a goal of enhancing the vehicle design 
and developing a capability for investigating 
autonomous and collaborative control technologies. 

Reference 4 documents the results of a wind tunnel 
test in which aerodynamic data was collected to provide 
a database to support the development of a dynamic 
simulation of the University of Florida MAV 
(UFMAV) concept. In that paper the flexible 
membrane wing was shown to significantly increase 
the stall angle of the vehicle without sacrificing L/D 
ratio. The vehicle was also determined to be statically 
stable in all axes. 

This paper describes the development of a dynamic 
simulation and flight control assessment based on the 
aerodynamic data described in reference 4. A 
control/guidance system design is also presented. The 
inner loop controller design uses measurement-based 
nonlinear dynamic inversion. The structure of the 
guidance system allows the vehicle to be integrated 
into an existing multiple vehicle collaborative control 
scheme. 151 


Table 1 - UFMAV geometric and mass properties. 


Empty Weight 

0.12 lbs 

Wing Area 

19.8 in 

Span 

6 in 

Mean Chord 

3.3 in 

Moments of Inertia: 


Ixx 

0.086 lb -111 

Iyy 

0.23 lb-in 

Izz 

0.21 lb-in 

Ixz 

0.037 lb-in 


Vehicle Description 

The University of Florida MAV (UFMAV) 
incorporates a high mounted flexible membrane wing 
and low mounted cruciform tail attached to a tapered 
fuselage with rectangular cross section (see figure 1). 
The fuselage is a truss-like design constructed of a 
carbon fiber/epoxy material covered with a thin 
transparent monofilm membrane. A more detailed 
description of the vehicle and its construction can be 
found in reference 3. Table 1 summarizes the pertinent 
geometric and mass properties of the vehicle. 

A unique aspect of the vehicle is its flexible 
membrane wing. The cambered wing structure is 
constructed of unidirectional carbon fiber prepreg 
laminate forming a leading edge spar and chordwise 
ribs or battens. A membrane material is bonded to the 
spar and batten. The wing membrane material is a 4 
mil thick flexible latex membrane. 

The maximum dimension (including length and 
wing span) of the vehicle is six inches. The wing area 
is approximately 19.8 square inches. The root chord is 
4.25 inches and the mean chord is 3.3 inches. The 
camber of the unloaded wing is approximately 6.5 
percent of the root chord with the maximum camber 
occurring at approximately 30 percent chord and is 
uniform across the span. The wing is mounted at an 
incidence of approximately nine degrees with wing 
incidence defined as the angle between the root chord 
line and the longitudinal axis of the fuselage. 

Control is accomplished using two independently 
controlled elevons that are actuated symmetrically and 
antisymmetrically using small rotary servos. A small 
gas engine normally provides propulsion with a three 
inch diameter propeller with a pitch of 1.25. However, 
an electric motor was used during wind tunnel tests to 
more accurately control propeller rpm and is used in 
the simulation model as well. 
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Simulation Model 

The simulation model is based on the aircraft 
equations of motion presented in Reference 6. The 
equations of motion were coded using 
Matlab/Simulink. [71 The structure of the simulation is 
depicted in Figure 2. The simulation is structured 
using subsystems representing actuator dynamics, 
equations of motion (EOMs), and sensor dynamics. A 
more detailed block diagram appears in the Appendix. 

The actuators subsystem currently consists of first 
order actuator transfer functions and limiters that 
bound the permissible range of symmetric ft S ym and 

antisymmetric 8 asy control surface deflections (±25 

degrees and ±20 degrees, respectively) and commanded 
motor voltage (0 - 20 volts). The sensor subsystem 
currently contains no dynamics but will permit sensor 
models to be added at a later time. 

The equations of motion include the longitudinal 
and lateral-directions equations of motion, models for 
thrust and aerodynamic forces and moments, and a 
standard atmosphere model (see Figure 3 and the 
Appendix). Each of the major components of the 
EOMs subsystem will be described subsequently. 
Equations of Motion 

The equations of motion are implemented in two 
major subsystems representing the vehicle dynamics in 
the longitudinal and lateral-directional axes. There is 
coupling between these two subsystems due to inertial 
and gravitational coupling. There are also several 
quantities that are used to determine the aerodynamic 
forces and moments (e.g., body rates, angle of attack, 
sideslip angle, speed). These quantities are fed back to 
the aero model as necessary. 



Figure 2 - UFMAV simulation structure. 



Figure 3 - EOMs subsystem structure. 


The longitudinal and lateral-directional subsystems 
consist of additional subsystems that systematically 
build up the equations of motion as derived in 
reference 6. The equations represent the six degree-of- 
freedom motion of a rigid aircraft relative to a flat, 
non-rotating earth. The atmosphere is represented using 
the 1976 Standard Atmosphere model. 161 

UFMAV Aero Model 

The aerodynamic model was obtained primarily from 
wind tunnel data collected in the NASA Langley Basic 
Aerodynamics Research Tunnel (BART). 1 Linear 
regression analysis was used to generate functions that 
approximate the dependence of the forces and moments 
on angle of attack, sideslip angle, and propeller rpm. 
The functions are in the form of Taylor series. 

The regression analysis was performed using wind 
tunnel data that consists of the aerodynamic force and 
moment coefficients at various combinations of angle 
of attack, sideslip angle, control surface deflection, 
dynamic pressure, and motor rpm. The range of 
variation for these parameters correspond to the region 
over which the aerodynamics are linear. 141 The main 
implication of this simplification is that the angle of 
attack is limited to values below 20 degrees and 
sideslip to values between -5 and 5 degrees. Cross 
terms between angle of attack, control deflection, and 
motor rpm are used to account for the dependence on 
propeller slipstream effects and the effect angle of 
attack has on control effectiveness. The values of the 
coefficients are shown in the tables in the Appendix. 

Note that there are three sets of coefficients for lift, 
drag, and pitching moment. Each set corresponds to a 
different dynamic pressure. The differences are 
attributable to Reynolds number effects. Interpolation 
is used in the simulation to determine the coefficient 
values at any given dynamic pressure between 1.0 and 
2.0 psf. Lack of sufficient lateral-direction force and 
moment data and higher levels of uncertainty for these 
quantities made it impossible to isolate Reynolds 
number effects for side force, rolling and yawing 
moment coefficients. As a result, the values for the 
lateral-directional coefficients represent an average over 
dynamic pressure. 

Additional terms were added to the Taylor series in 
an ad hoc manner to account for dependence on angular 
rates (i.e., dynamic derivatives). Terms associated with 
the angle of attack and pitch rates were added for lift 
(Ci^,Ci.) and pitching moment (Cm q ,CM & )- 

Terms associated with roll and yaw rates were added for 
the side force (Cy^ , Cy r ), rolling moment (C;^ , 

Ci r ) and yawing moment (C n ^,C„ r ). The 
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coefficients for most of these terms were computed 
using PMARC. [81 Two exceptions were the dynamic 
derivatives associated with rate of change of angle of 
attack (Ct ,Cm ) which were chosen based on 

“typical" values published in reference 9, page 19. The 
values for all the dynamic derivatives are shown in the 
tables in the Appendix. 

The expression for total lift force coefficient is 
shown in equation (1) as an example of the Taylor 
series expansion. 


Ci = K+CiU) + K +c Ii<r r) 

Ct. + C/- 7 T 

Osym osym* I y 
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in 


+Cr 


( \ 
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y2V y 
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q+C L 
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The thrust dependent cross terms C/, account for 

the fact that the effects of thrust are coupled with 
angle of attack and control surface deflection through 
prop stream effects. The expressions for drag, side 
force, and pitching, rolling, and yawing moments are 
similar in structure but differ in the particular 
coefficients associated with coupling. 

UFMAV Propulsion Model 

The propulsion model was obtained from wind 
tunnel data collected during the BART test. 141 Motor 
thrust was approximated by subtracting the prop-off 
axial force from the prop-on axial force. Regression 
analysis was used to generate generalized Taylor series 
functions that approximate the dependence of motor 
thrust on angle of attack, dynamic pressure, and 
voltage commands. The propulsion model consists of 
two parts: a motor model that characterizes the 
relationship between motor voltage command and 
propeller rpm, and a thrust model that characterizes the 
relationship between propeller rpm and thrust 
coefficient. This implementation separates the effect of 
propeller loading on motor rpm from the thrust 
produced at a given rpm. 

The regression analysis was performed in an ad hoc 
manner to identify a combination of parameters that 
provide a reasonable approximation to experimental 
data. The function approximating the relationship 
between motor voltage and motor rpm is 

KRPM= C Mq + C. M q + C Ma a +C M or 

y a 

2 ^ 

v motor ^ v motor 

v motor v , 

motor 


where KRPM represents rpm/1000. Note that the 
behavior is essentially quadratic in motor voltage 
( v mo t or ) with a variable offset which is determined 
by propeller loading effects expressed through a 
dynamic pressure ( q ) dependent term and the angle 
of attack dependent terms. 

The function relating thrust to motor rpm is 


C. T = (Cy, +Ct 2 KPRM ) 


^krpm! 


(3) 


where Cj is the thrust coefficient. No attempt was 
made to explain the structure of this equation on a 
physical basis. Tables of the propulsion model 
coefficients are presented in the Appendix. 


Analysis 


The simulation model of the UFMAV was used to 
perform a number of analyses to assess the stability 
and control properties of the vehicle. These analyses do 
not, however, constitute a validation or verification of 
the simulation model since there are no static or 
dynamic data available for the actual aircraft in flight. 

First a trim comparison is made for the vehicle in 
straight level flight at several dynamic pressures. The 
dynamic pressures (1.0, 1.6, and 2.0 psf) correspond to 
conditions at which experimental data are available. 
These data were obtained during the wind tunnel test in 
BART 141 and are representative of typical flight speeds 
of the UFMAV. 

The results of three longitudinal trim studies are 
shown in table 2. The experimental trim results were 
obtained by achieving trim in the BART tunnel. This 
was accomplished by first setting the tunnel speed 
corresponding to the desired dynamic pressure and then 
varying the vehicle angle of attack, symmetric elevon 
deflection, and motor voltage (i.e., propeller rpm) until 
the lift was approximately equal to the gross vehicle 
weight and the pitching moment and total axial force 
were both approximately zero. 

The simplified analytical trim was determined using 
the method described in reference 10. Equation (4) is 
the matrix equation that was solved to determine trim 
angle of attack and symmetric elevon. 


C, 
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Mr, 
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sym 


M S 


sym 


[ OC trim 
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w c 
— 7 “ '-'Ln 

■ = 4 

qS (, \ 
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l 
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The lift curve slope, moment curve slope, and lift and 
moment control sensitivities were obtained from the 
experimental data for the corresponding dynamic 
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Table 2 - Experimental, analytical, and simulation 
based longitudinal trim. 


Dynamic 
Pressure 
fesf) 

Propeller 

RPM 

Angle of 
Attack 
(deg) 

Symmetric 

Elevon 

(deg) 

Experimental Trim 

1.0 

18,900 

10.4 

-6.5 

1.6 

20,600 

5.4 

-3.5 

2.0 

21,900 

4.0 

-2.5 

Analytical Trim (Simplified) 

1.0 

- 

11.2 

-5.6 

1.6 

- 

5.4 

-2.5 

2.0 

- 

3.5 

-1.9 

Computed Trim (UFMAV) 

1.0 

19,600 

11.1 

-6.8 

1.6 

21,200 

5.6 

-4.7 

2.0 

22,000 

3.5 

-1.9 


Table 3 - Computed lateral-directional trim. 


Dynamic 

Sideslip 

Bank 

Antisymmetric 

Pressure 

Angle 

Angle 

Elevon 

(psf) 

(deg) 

(deg) 

(deg) 

Computed Trim (UFMAV) 

1.0 

-0.051 

-0.97 

-0.59 

1.6 

0.028 

-1.6 

-0.54 

2.0 

0.070 

-2.1 

-0.51 


pressure. 141 The data was assumed to correspond to a 
propeller rpm near trim. 

The computed trim was obtained by using the 
UFMAV simulation model and a constrainted 
optimization routine to achieve level trim at a specified 
dynamic pressure. Comparison of the three trim 
analyses shows very good agreement for angle of 
attack, symmetric elevon deflection, and propeller rpm. 
This implies that the longitudinal aerodynamic forces 
and moments are well approximated in the simulation. 

A straight and level trim analysis using the 
simulation model was also performed to determine the 
lateral-directional quantities: sideslip, bank, and yaw 
angles. Table 3 shows the results of this analysis. 
Note that the UFMAV achieves lateral- directional 
control via antisymmetric elevon and dihedral 
coupling. It does not have two independent lateral- 
directional controls (such as rudder and aileron) and 
cannot be trimmed at zero bank angle (or zero sideslip 
angle) as is typical. The results indicate that though 
the vehicle does have significant asymmetries, all the 
trim values are small and within the range of values at 
which the aerodynamic data was obtained and are 
qualitatively consistent with the vehicle in flight. 


Table 4 - longitudinal modes. 


Dynamic 

Pressure 

(psf) 

Short Period Mode 

Phugoid Mode 

damping 

ratio 

freq. 

(rad/sec) 

damping 

ratio 

freq. 

(rad/sec) 

1.0 

0.13 

23.3 

0.44 

0.85 

1.6 

0.12 

30.2 

0.35 

0.65 

2.0 

0.12 

32.6 

-0.56 

0.67 


Table 5 - lateral-directional modes. 


Dynamic 

Pressure 

(psf) 

Spiral 

Mode 

Roll 

Mode 

Dutch Roll Mode 

eigenvalue 

eigenvalue 

damping 

ratio 

freq. 

(rad/sec) 

1.0 

-1.04 

-27.7 

0.094 

21.1 

1.6 

-1.04 

-37.3 

0.065 

24.2 

2.0 

-1.02 

-42.8 

0.050 

25.9 


The simulation was also linearized about the above 
trim conditions to assess the dynamic stability of the 
vehicle. Table 4 summarizes the frequency and 
damping of the linearized longitudinal modes. Note 
that the short period mode is stable for all three 
dynamic pressures but lightly damped. Its frequency 
increases with increasing dynamic pressure but the 
damping is essentially constant. The damping of the 
phugoid mode varies significantly and is unstable at 
the higher dynamic pressure. 

Table 5 summarizes the eigenvalues or frequency 
and damping of the linearized lateral-directional modes. 
Note that all the modes are stable and that the dutch 
roll mode is lightly damped. This is qualitatively 
consistent with behavior of the vehicle in flight. Note 
that the spiral mode is relatively unaffected by changes 
in dynamic pressure but that the magnitudes of both 
the roll and dutch roll modes increase with increasing 
dynamic pressure. 

Linearized models used to perform this analysis can 
be found in the Appendix. 

Control Design 

A preliminary guidance/control system has been 
developed to enable investigations of autonomous and 
collaborative control issues. The controller is 
composed of two main parts: an inner-loop 

measurement-based nonlinear dynamic inversion 
controller for control of angular rates and an outer-loop 
navigation command follower for control of wind-axis 
angles. 111 ' 1 " 1 An overview of the control system is 
given in figure 4. The control system inputs are 
commanded flight -path angle y, wind-axis heading 
angle %, and total speed V t . These inputs were chosen 
to allow the vehicle to be readily integrated into an 
existing multiple vehicle collaborative control 
scheme. 151 Controller outputs are commanded symmetric 
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Figure 4 - structure of UFMAV control system. 

and antisymmetric elevon deflection. A separate 
proportional-integral error loop is used to generate 
motor voltage commands to control total velocity. 
For this preliminary study, the feedback 
measurements are assumed to be known perfectly. 
The two main parts of the control system are 
discussed in more detail in the following. 


Measurement-based Nonlinear Dynamic Inversion 
Given desired values of roll acceleration p, pitch 
acceleration q , and yaw acceleration r, the inner- 
loop controller generates symmetric and 
antisymmetric elevon commands to achieve the 
desired angular accelerations. The inner-loop 
controller is based on a modified nonlinear dynamic 
inversion approach developed in reference 1 1 . This 
approach does not require a model of the baseline 
vehicle (i.e. no stability derivatives), but does require 
a model of the vehicle's control effector derivatives 
and feedback of body-axis angular accelerations and 
control effector positions. Since this approach uses 
acceleration measurements in lieu of a complete on- 
board vehicle model, this approach is less sensitive to 
vehicle model errors and can adapt to vehicle failures 
and/or damage. An overview of the approach from 
reference 1 1 is given in the following. 

Given the vehicle equations of motion 

x=F(x,8)=f(x)+g(x,8 ) 

T (5) 

y = [p q r] =h(x) 


where x is the vehicle state vector, § is the vehicle 
control vector, and y is the vector of control 
variables: roll rate p, pitch rate q, and yaw rate r. A 
Taylor series expansion of (5) yields the following 
first-order approximation to F(x, 8) in the 
neighborhood of [.Tq ,§o] 


F(x,8)=f(x 0 )+g(x 0 8 0 )+ 


(f(x)+g(x,8)) 


Ts( s(x - 5> > 


(.r-.ro) 


x=.r 0 ,5=§o 

(5 -So) 


(6) 


x=Xq 8 


Letting x 0 and 8 0 denote a previous state and control 
from the recent past and defining 


•^° = — ( /(r)+g(r,8) } | 


5 ° = ^:(g(*>S))| 


x=x 0 $=S 0 
x=x 0 ,8=8q 


(7) 


F( x, 8) can be written as 

F(x, 8) = x = xq + Aq(x- a'o ) + 5o A8 (8) 

in the neighborhood of X—Xp, 8=8o where 
8 = 8 0 + A8. 

At this point, this development differs from 
reference 1 1 in that the number of controls is less than 
the number of controlled variables and so the desired 
responses cannot be completely achieved. A control law 
is obtained by minimizing 

J = i}’d~y) T Q(Vd-y) (9) 

where 


y = 



h x (xo +A 0 (x-xo) + B 0 A8) 


(10) 


and Q is a positive-definite diagonal weighting matrix 
used to emphasize desired system responses. This 
yields 


A5 = 


(hxBo) 1 Qh x B(, 


( More • 


(Vrf-M o) 


(ID 


With a sufficiently fast update rate x tends to x 0 and 
equation (11) becomes 


A8 = 


i-l 


{h.xBof Qh x B 0 ( h x B 0 ) T Q(y d -h x xo ) 


( 12 ) 


where 8 = 8 0 + A8. The vehicle's control derivatives 
B 0 are generated from the nonlinear aerodynamic 
control coefficients using a central difference 
approximation. 


Navigation Command Follower 

Given desired values of flight path angle y, wind- 
axis heading %, and total velocity V t the navigation 
command follower generates required roll rate, pitch 
rate, and yaw rate acceleration commands for the 
inner-loop controller. 

The desired dynamics for the outer-loop were chosen 
to be 
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id =co t (Yc -y) 

id=<* x (Xc-x) 


where the subscript d denotes the desired value and the 
subscript c denotes commanded input values. The 
bandwidths co x and co 7 were chosen to be 
approximately a decade below the bandwidths of the 
desired inner-loop dynamics and therefore were 
chosen to be 2 rad/sec. 

Using the wind-axis point mass equations of motion 
and assuming sideslip angle and sideforce are small and 
that V t and cosfy) are non-zero, commanded wind- 
axis bank angle g c can be determined as a function of 

v t>v d and x rf 1131 


tan \i c = 


V,XdC osy 

V,y d + gcosy 


(14) 


The desired dynamics for wind-axis bank angle was 
chosen to be 

M-rf=<M|4 c -|4) (15) 

where co^ were chosen to be 4 rad/sec. 

The wind-axis angular rates (I, y, and % are 

transformed to commanded body-axis rates (assuming 
sideslip angle is zero) using 


Pc 


cosa 0 -sina 

1 0 -siny 


dc 

> — 

0 1 0 

0 cosg sing cosy 

Y 

r c_ 


sina 0 cosa 

0 -sing cosgcosy 

X 


where a is angle of attack. The desired closed-loop 
dynamics y for the inner-loop were chosen to be 

Pd= a p(Pc~P) 

id = (Oq{q c -q) (17) 

id =(O r ( r c ~ r ) 

where the subscript d denotes the desired value and the 
subscript c denotes commanded values determined by 
the outer-loop control law. The inner-loop 
bandwidths co p , co q , co r were chosen to be 20, 15 and 
20 rad/sec, respectively, consistent with the open- 
loop bandwidth. 

Figure 6 shows time responses for inner-loop p 0 , q c , 
and r c angular rate commands (i.e. no outer-loop 
controller). The initial condition for these time 
responses is straight and level flight at V t = 37 feet/sec. 
Reference signals were generated for comparison with 
the achieved responses from equations (13), (15), and 
(17) using the specified bandwidths. The commands 
and reference signals are shown respectively as dashed 
and dotted lines in the upper two plots of figures 6a 


and 6b. A stability-axis roll rate doublet was 
commanded (50 deg/sec from 1 to 2 seconds and -50 
deg/sec from 2 to 3 seconds) with pitch rate and 
stability-axis yaw rate commanded to zero. The effect 
of choice of control variable weighting is demonstrated 
in these figures. Figure 6a shows responses for a 
control variable weighting of Q = diag([roll 
acceleration error weighting, pitch acceleration error 
weighting, yaw acceleration error weighting]) = 
diag([ 1,5,2]). As can be seen, the stability-axis roll rate 



Time (sec) 


Figure 6a - inner-loop p c , q c , and r c angular rate 
commands, Q = diag([ 1,5,2]). 



Time (sec) 


Figure 6b - inner-loop p c , q c , and r c angular rate 
commands, Q diag([ 1 ,5, 1 0]). 
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revolts ™ — ***~ 


Is” 


Time (sec) 

Figure 7 - Figure 3 - outer-loop Y c , % c > an d 
V( commands, Q = diag([ 1,5,2]). 


reference signal is more closely followed than the 
stability-axis yaw rate reference signal with low dutch 
roll damping as illustrated by the oscillatory sideslip 
response. Figure 6b shows responses for Q = 
diag([ 1,5,10]). This results in the stability-axis yaw- 
rate reference signal being more closely followed 
than in figure 6a and a better damped sideslip 
response. 

Figure 7 shows time responses for outer-loop Y c , 
X , and Vt commands. The control variable 

weighting was Q = diag([l,5,2]). The commanded and 
reference signal values are shown respectively as dashed 
and dotted lines in the top three plots. As can be seen, 
the vehicle closely follows the reference signal with 
reasonable control activity. 

These preliminary results demonstrate that this is a 
viable approach for control of systems where the 
number of controls is less than the number of control 
variables, such as, MAV's. Future efforts will focus on 
improvements to this approach, robustness analysis, 
and use of this method as part of a multiple vehicle 
collaborative control scheme. 

Concluding Remarks 

A dynamic simulation model of an aeroelastic fixed 
wing micro aerial vehicle has been developed that is 
suitable for a wide variety of uses including control 
system design, navigation and guidance algorithm 
development, and their assessment. The simulation is 
based on a vehicle concept developed at the University 
of Florida and wind-tunnel data collected in the NASA 
Langley Basic Aerodynamics Research Tunnel. 
Regression analysis was used to obtain a generalized 
Taylor series aerodynamic model. 


The simulation was used to assess vehicle trim and 
basic stability and control properties. The analysis 
indicates that the vehicle has acceptable stability 
properties and good controllability. 

A control system was designed using a 
measurement-based nonlinear dynamic inversion 
approach. The method was extended to accommodate 
application to systems with fewer controls than 
controlled variables as is the case for the subject 
vehicle. A guidance loop was also designed to allow 
the simulation model to be integrated into an existing 
multiple vehicle collaborative framework. 

Assessment of the control and guidance systems 
using the simulation demonstrated satisfactory 
performance. Additional research is underway to 
improve the dynamic response, investigate performance 
robustness, and explore implementation issues. 
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Appendix 

The appendix contains a block diagram of the basic 
structure of the simulation model, tables of the 
aerodynamic force and moment coefficients and 
linearized models of the UFMAV. 


Simulation Block Diagram 
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Aerodynamic Force/Moment Coefficients 
Longitudinal Axes 


Lift 

Coefficient 

mm 

m 

q = 2.0 
(psf) 

Cl 0 

0.36094 

0.341 

0.344 

C L a 

3.1365 

2.95 

2.94 

C L * 

°sym 

0.65369 

0.590 

0.590 

Cl l 

0.0043367 
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0.00169 

C LuT 

0.035346 
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0.0228 

^§syniT 

0.012241 
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CL q 

4.0 

4.0 

4.0 

C L d 

1.0 

1.0 
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Drag 

Coefficient 

q = 1.0 
(psf) 

q = 1.6 
(psf) 
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C D a 
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°sym 
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°sym 
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C Di 
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2.99e-4 

8.24e-5 
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C °ah 

ao sym 

0.648 

0.731 

0.659 


Pitch Moment 
Coefficient 

q= 1.0 
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q= 1.6 
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q = 2.0 
(psf) 

Cm o 

0.0797 

0.0321 

0.0304 

C »>a 
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-0.702 

-0.741 

Cm* 

°sym 

-0.350 

-0.384 

-0.355 

Cm l 

2.61e-4 

1.26e-4 

-4.86e-5 

m aS S ym 
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0.0884 
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C »>aT 
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CjYl s: 

° sytm 

-5.5 le-3 

-4.52e-3 

-2.52e-3 

Cm q 

-0.5 

-0.5 

-0.5 

C m • 

m a 

-3.0 

-3.0 

-3.0 


Lateral-Directional Axes 


Coefficient 

Sideforce 

(Y) 

Rolling 

Moment 

(1) 

Yawing 

Moment 

(n) 

C «o 

-1 .9 le-2 

-4.72e-3 

5.73e-3 


-8.55e-l 

-2.74e-l 

-3.25e-l 

C(,) 5 asy 

-5.72e-2 

6.00e-2 

3.00e-2 

C O r 

1.83e-3 

2.70e-4 

-2.46e-4 

C( *^asyT 

-3.57e-3 

1.25e-3 

1.82e-3 

< '(*)(37' 

-2.15e-2 

-4.09e-3 

-3.33e-3 

C{ *\bodv 

-0.35 

-0.65 

-0.015 

C (*\ody 

1.3 

-0.32 

-0.5 


Motor and Thrust Models 


Motor 

Coefficient 

Value 

Cm o 

2.17 

c M- q 

0.0793 

v motor 

-0.0620 

Cm 2 

v motor 
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c M a 

-0.0264 

Cm ? 

a 

-0.179 


Thrust 

Coefficient 

Value 

Ct\ 

-0.00140 

Ct 2 
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Linear Models 
Longitudinal Axis 

T 

input vector: u = [O^ v motor \ 

T 

state vector: x = [V( « q 0] 

q = 1.0 psf 


-1.17 

-11.3 

0 

-32.2 

-0.0602 

-4.02 

0.982 

2.99e-5 

-5.14 

-545. 

-1.80 

0 

0 

0 

1.0 
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-0.121 

0.00136 
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1.61 



B= 
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0.624 



B= 

-0.852 

-0.0193 
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-0.179 



-327. 

-0.764 




0 

0 




0 

0 







.6 psf 




q = 1.6 psf 





" -0.648 

-6.55 

0 

-32.2 


-1.58 

0.0954 

-0.985 

0.873 


-0.0355 -4.67 

0.982 3.40e-5 


-2020 

-28.5 

-17.8 

0 


-4.43 

-915. 

-2.28 

0 

A = 

261 

-5.27 

-11.4 

0 


0 

0 

1.0 

0 


0 

1.0 

0.0987 

0 



"-3.76 

1.23 




-0.160 

0.00148 



B= 

-0.896 

-0.00814 



B= 

632 

0.893 




-573. 

0.238 



277 

-0.270 




0 

0 




0 

0 



q = 2.0 psf 


q = 2.0 psf 


f -0.670 

-4.5 

0 

-32.2 


—1 .79 

0.0587 

-0.987 

0.782 

-0.0247 -4.87 

0.982 3. 14^-5 


-2562 

-31.2 

-19.9 

0 

-2.38 

-1067. 

-2.55 

0 

A = 
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-5.89 

-12.8 

0 

0 

0 

1.0 

0 


0 

1.0 

0.0617 

0 


’-3. 12 

4.16 



-0.183 

0.00148 


B= 

-0.912 
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-0.00428 

0.161 


B= 
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1.07 

-0.334 



0 

0 



0 

0 



Lateral-Directional Axis 

T 

input vector: u =[S asy v motor] 

T 

state vector: x = [ P P r §\ 

q = 1.0 psf 


-1.22 
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-0.971 

1.09 
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-14.1 

0 
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-4.17 

-9.01 

0 

0 

1.0 

0.196 
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